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Abstract 

We present an analytically solvable model of P colinear, two-dimensional 
quantum dots, each containing two electrons. Inter-dot coupling via the 
electron-electron interaction gives rise to sets of entangled ground states. 
These ground states have crystal-like inter-plane correlations and arise dis- 
continously with increasing magnetic field. Their ranges and stabilities are 
found to depend on dot size ratios, and to increase with P. 
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Recent experimental and theoretical interest in quantum dot systems has opened up 
the fascinating area of highly-correlated, few-body quantum phenomena in the traditionally 
large-A^ field of semiconductor physics. Complex ground state behavior as a function of 
magnetic field has been predicted for single two-dimensional (2D) quantum dots containing 
N electrons [|^. Even for as few as = 2 electrons per 2D dot, "magic number" ground state 
transitions are predicted as a function of magnetic field as a result of the competition between 
(single-electron) confinement energy and (many-body) electron-electron interactions P|. 
Remarkably, such transitions for = 2 have recently been observed experimentally [Q. 
Arrays of coupled dots @] [H have been attracting increasing attention, partly because 
of the possibility of apphcation as ultra-small logic gates [|1^ [0. Adjacent dots can 



be coupled by "optical wiring", i.e. coupled by the two-body electron-electron interaction 
between electrons in adjacent dots which can be non-zero even in the absence of a single- 
body tunneling term [|l^ . Such few-electron coupled 2D systems are also interesting in that 



they can represent the small- A^ analogs of coupled, parallel 2D electron gases ||T3|. Given 
the complex ground state behavior in a single 2D dot as a function of 5-field, the coupled 
2D dot system offers the interesting possibility of competition between electron correlations 
on the same dot (i.e. intra-dot interactions Vintra) and correlations between adjacent dots 
(i.e. inter-dot interactions Vinter)- 

This paper predicts the existence of entangled ground states with crystal-like inter-dot 
correlations in multiple quantum dot systems. These ground states occur discontinously 
with increasing magnetic field, being interdispersed with states having negligible inter-dot 
correlations. The ranges and stabilities of the crystal-like states depend on dot size ratios 
and increase with the number of dots P. Our model is solved analytically and consists of 
P colinear 2D quantum dots, each containing two electrons but not necessarily identical 
in size. The model considers a sufficient number of electrons as to contain both inter and 
intra-dot electron-electron interactions, and yet still admit analytic solutions. These analytic 
solutions implicitly include mixing with all Landau levels. 

Figure 1 shows our system for a pair of dots (P = 2). The dots are arranged vertically 



with separation s. This vertical geometry is of specific experimental interest given the 
possibility of fabrication via etching of a multiple quantum well structure (see Ref. for 
P = 2). Following several single-dot studies we model each of the P dots by 2D (xy 
plane) parabolic potentials with a perpendicular magnetic field B (z direction) of sufficient 
strength to spin-polarize the electrons. Electrostatic confinement in the z-direction is taken 
to be sufficiently strong that the electrons are frozen in the lowest z sub-band. We take 
the electron-electron interaction potential to be of inverse- square form, ^ where /? is a 
positive parameter; this interaction for a single layer gives quantitatively similar results to 



the bare Coulomb interaction [T^ ||T6[. Inter-dot tunneling is assumed to be negligible. The 



Hamiltonian for the system of P dots is (with a symmetric gauge) H = Hq + V with 

a=l,2,...P;i=l,2 ^ ^ 

Pa,i and Ya^i are the momentum and position of electron i in dot a. Each electron has effective 
mass m* and z-component of angular momentum Lq, j. The cyclotron frequency is uo^ and 

uj^[B) = uJq^ + ojI/A. The electrostatic confining potential oj^a is in general different for 

_ 1 

each dot (N.B. the dot 'size' ~ c^oa )• The dominant inter-dot coupling is due to interactions 
between electrons on adjacent dots; we therefore take V = Vintra + Vinter where 

Vintra = E Z 12 



a=l,2,...P 

and 



Vinter | Z ~ 12 , 2 ' (^) 

a=2,3,...P i,j 

For Vinter ~^ 0, the exact eigenstates of H are products of P single-dot states [|T^. For finite 



Vinter, HO cxact analytic solutions are known. Our approach is to Taylor expand Vinter under 
the assumption that \ra,i — r(a_i)jp < s^. Having solved the resulting problem, we can 
check from the analytic expressions for (e.g.) ^ that any particular set of dot parameters 
is consistent with this assumption. Hence 

Vnter = ^ E E E (-1)^^^^^^^^^^)' (4) 



Exact analytic solution of H is now possible including terms of order (3r'^/s^ in Vinter', (ana- 
lytic) perturbation theory is then employed for the /5r^/ terms. 

We now discuss the method explicitly for P = 2. Employing an orthogonal transforma- 
tion, with coefficients depending on the relative dot sizes uji and 002, allows exact solution of 
H including terms of order /3r'^/s'^ in Vinter- The only non-trivial contributions to H to this 
order are those which depend on the relative positions of the electrons within each dot. The 
corresponding eigenstates contain quantum numbers rria (the relative angular momentum 
between the two electrons in dot a). The ground and low-lying states of the P = 2 system 
have all other quantum numbers zero; these states can be labelled \mi,m2) (signifying the 
direct product of |mi) and 1^2)) and are used as a basis for the Pr^/s^ perturbation. Figure 
2(a) shows the ground state transitions as a function of -B-field in the limit s ^ 00 (i.e. 
Vinter 0) for the case of equal dot sizes (tui = 002). As a demonstration of the (lack of) 
inter-dot correlation, the insets show the charge density in a given dot with the electrons in 
the other dot 'fixed' opposite one another at the crosses. The plots are angularly symmetric 
(i.e. negligible inter-dot correlation). The radial localization increases as the ground state 
rria increases (i.e. as B ^ 00, each dot approaches its own classical limit). Figure 2(b) 
shows the corresponding diagram for finite s. The main feature is that new entangled |T^ 
states arise as a result of finite V^ter (i-e. including terms to order f3r^/s^). The density 
plots show the strong inter-dot correlation which characterises the new entangled ground 
states. 

These entangled states can be easily understood for equal dot sizes. First note that in the 
exactly solved (i.e. order Pr^/s^) system the states \ma,mi,) and \mh,ma) are degenerate. 
Secondly, the perturbation (i.e. order f3r'^/s^ terms) only mixes \ma,mb) and \mc,md) if 
"^a — = ~ = or ±2. The mixed (i.e. entangled) states are therefore (|m -|- 
2, m) ± \m, m + 2)) which we write as \m + 2, m)+ and \m + 2, m)_. It is remarkable that as 
B increases the ground state switches back and forth between 'pure' states (mi = 1112) and 
entangled states (N.B. the |m,m) states are now not strictly 'pure' because of a very small 
non-degenerate perturbation mixing). For large B {B > lOT), the entangled states prevent 
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the pure states from becoming ground states. As S — > oo, the classical hmit is reached of 
four point charges situated at the corners of a square when projected onto the xy-p\ane. An 
alternative view of the formation of these entangled states is as a resonance phenomenon; 
for example, the state |3, 1) can be thought of as continually exchanging energy with |1, 3) 
via virtual photons (i.e. via the electromagnetic field representing the electron-electron 
interaction). This is equivalent to the resonant Forster process which is well-known as an 
energy transfer mechanism for biological molecules and proteins. In excitonic language, 
the states |3, 1) and |1,3) correspond to the |1, 1) vacuum plus an electron-hole excitation 
(exciton) of angular momentum 2 on dots 1 and 2 respectively; the formation of the entangled 
state |3, 1)_ represents the resonance between these two adjacent excitons. 

Figures 3(i)a and 3(i)b show in detail the region near the |1,1) and |3,3) crossover 
discussed above. Figures 3(ii)a and 3(ii)b show the corresponding region for the case of 
unequal dot sizes {cui ^ lo-i). In this case |ma,m6) and |m6,ma) are no longer degenerate. 
However by keeping \lo\ — lo-2\ « cui, we can arrange that the two states are energetically 
much closer to one another than to any other state with which they can mix. Under these 
conditions we can again solve the problem analytically; the entangled states are (p|m -|- 
2,m) ± q\m,m + 2)) where p,q are unequal and depend on the energies of the component 
states. Figure 3(ii)b shows that the stability of the entangled ground states, i.e. the gap 
between the entangled ground states and the other competing states, is greater than the 
corresponding gap in Fig. 3(i)b. The gap scales approximately as \Je^ -\- {/S.Ey where e 
is the gap in the equal dot system (Fig.3(i)b) and Ai? is the energy difference between 
component states (i.e. the separation of the dashed lines in Fig. 3(ii)a). The results of Figs. 
3(i) and 3(ii) are analogous to linear and quadratic Stark effects; in both Figs. 3(i) and 3(ii) 
we essentially have two-level systems with the former being degenerate while the latter is 
non-degenerate. 

Figures 3(iii)a and 3(iii)b consider the case of P = 3 equal-size dots. Exact analytic 
treatment to order /3r^/s^ proceeds as above. The ground and low-lying states are now 
characterised by three non-zero quantum numbers describing the relative angular momentum 
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between electrons on each of the three dots; we label these states |mi, 777,2,7713). Just as for 
P = 2 with cui = CU2 to this order (i.e. /Sr^/s"^) we find degeneracies: \ma,mb,mc) is 
degenerate with \mc,mb,ma). In the limit of large s all other states formed by permuting 
majmbjUic would also be degenerate. Having exactly solved for P = 3 to order /3r^/s^, we 
again turn to perturbation theory to treat the terms of order /3r^/ (which we will denote 
by h) . The non-zero matrix elements of h in the nearly degenerate subspace are found to be 
of the form (777a + X, 7775 + Y, rric + Z\h\ma, mj,, where {X, Y, Z} is any permutation of 
{—2, 0, 2}. For example, |3, 1, 1) mixes with |1, 3, 1) and |1, 1, 3) yielding three new mixed 
(i.e. entangled) states. These entangled states then compete with each other to become 
ground states for finite s (see Fig. 3(iii)b). The entangled states are the P — 3 dot analogs 
of the inter-dot correlated, crystal-hke states for P — 2. 

As the number of dots P increases, the gain in energy of the entangled (i.e. inter-dot 
correlated) states at finite s actually increases as compared to the s — > 00 limit. This increase 
is non-linear with P. Note that we can think of the earlier cases of P = 2 and P = 3 
as a 'diatomic molecule' and 'triatomic molecule' respectively. Consider the 'polyatomic 
molecule' with P identical 'atoms' (dots) all in the 777 = 1 state. The state m — 3, for 
example, can be created on any one of these dots yielding P degenerate 'tight-binding' 
combinations or 'molecular states' in the limit s — > 00. The degeneracy of these states is 
broken for finite s by the coupling between adjacent dots. In contrast to an actual tight- 
binding molecule where it is the single-body tunneling term that couples atoms, here the 
coupling is via the two-body interaction Vinter- The states represent 'Frenkel excitons'; the 
exciton (i.e. m — 3 state) on a given dot can transfer its energy resonantly to all members 
of the chain. For large P, the analysis is simplified considerably by introducing periodic 
boundary conditions. This removes 'end-effects', introducing translational symmetry into 
the problem and allowing application of Bloch's theorem. Consider a system of P dots 
with a lattice structure; if a;, = a; for i = 1, . . . P then this is a 'monatomic crystal'. If 
alternating cj's differ, we have a 'diatomic crystal'. The entangled states are now travelling 
wave excitations and yield an exciton band structure. 
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In summary, we have presented details of entangled ground states arising in colinear, 
multiple dot structures. The tendency for formation of such states increases with the number 
of dots. 

This work was supported by an EPSRC Studentship (S.B.) and the Nuffield Foundation 
(N.F.J.) . 
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Figure Captions 

Figure 1: 

Schematic illustration of the P = 2 dot system. Each dot contains two electrons. 
Figure 2: 

(a) Energies (in Kelvin relative to an arbitrary zero) of low-lying eigenstates | mi, 7712) as a 
function of S for P = 2 dot system in limit s — > 00. Lowest curve at a given B corresponds 
to the ground state. From left to right the solid lines represent successive ground states 
|1,1), |3, 3) and |5,5); dashed lines correspond to the degenerate pairs (|3, 1), |1,3)) and 
(|5,3), |3,5)). Contour plots are ground state charge densities in a given dot with electrons 
in the other dot fixed at the crosses. 

(b) As (a), but with finite s. The degenerate pairs from (a) spht to form entangled states; 
the dashed hues to the left represent new states |3, 1)+ and |3, those to the right are 
|5,3)+ and |5,3)_. States |3, 1)_ & |5,3)_ become ground states; their charge densities are 
shown. 

Figure 3: 

Energies (in Kelvin relative to an arbitrary zero) of low-lying states as a function of B for 
three different systems, (i) P = 2 dots of equal sizes. Parts (a) and (b) are magnifications 
from Figs. 2(a) and 2(b) respectively, (ii) P — 2 dots of unequal size, (iii) P = 3 dots of 
equal size; in (a) each dashed line is triply degenerate while in (b) only the lowest two of 
the six entangled states are shown. 
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